arXiv: 1506.02661 vl [hep-ph] 8 Jun 2015 


CERN-PH-TH-2015-134 

STUPP-15-222 


Effective field theory approach to 6 —)■ sU^'\ B —)■ and B 

with third generation conplings 

Lorenzo Calibbi/’^ Andreas Crivellin,^ and Toshihiko Ota^ 

^ State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, 

Chinese Academy of Sciences, Beijing 100190, P. R. China 
^Service de Physique Theorique, Universite Libre de Bruxelles, C.P. 225, B-1050 Brussels, Belgium 
^ CERN Theory Division, CH-1211 Geneva 23, Switzerland 
‘^Department of Physics, Saitama University, Shimo-Okubo 255, 338-8570 Saitama-Sakura, Japan 

LHCb reported anomalies in _B —>■ K*pApT, Bs —>■ p~ and R{K) = B —>■ Kp'^p~/B —> 
Ke'^e~. Fnrthermore, BaBar, BELLE and LHCb found hints for the violation of lepton flavour 
universality violation in R(D^*^) = B ^ D^ *\v/B —^ D^*^lv. In this note we reexamine these 
decays and their correlations to B —>■ using gange invariant dim-6 operators. For the 

nnmerical analysis we focus on scenarios in which new physics couples, in the interaction eigenbasis, 
to third generation quarks and lepton only. We conclude that such a setup can explain the b —>■ 
sp'^p~ data simultaneously with R{D^*^) for small mixing angles in the lepton sector (of the order 
of 7r/16) and very small mixing angles in the quark sector (smaller than Vcb)- In these region of 
parameter space B —>■ K^*'^Tp and Bs —>■ rp can be order 10“®. Possible UV completions are briefly 
discussed. 


I. INTRODUCTION 

So far, the LHC completed the standard model (SM) 
of particle physics by discovering the last missing piece, 
the Higgs particle Furthermore, no significant 

direct evidence for physics beyond the SM has been 
found, i.e. no new particles were discovered. However, 
LHCb observed indirect ‘hints’ for new physics (NP) in 
B -T K*p'^p~, Bs —)• (f>p'^p~ and R{K) = Br(B -T 
Kp~^P~)/By{B —>■ ATe+e”). Furthermore, BaBar and 
also very recently BELLE and LHCb reported lepton 
flavour universality violation in B —)■ These ob¬ 

servations can be used as a guideline in the exploration 
of possible physics beyond the SM. 

In more detail, the current experimental situation is as 
follows: LHCb reported deviations from the SM predic¬ 
tions [3] in B —>• K*p^p~ m |S] (mainly in an angular 
observable called Bg 0) with a significance of 2-3 a de¬ 
pending on the assumptions of hadronic uncertainties 0- 
m. Also in the decay Bg —>■ T - [lU] LHCb uncovered 

differences compared to the SM prediction based on lat¬ 
tice QCD mm and light-cone sumrules [13] of 3.1 cr 
mm Eurthermore, LHCb m found indications for the 
violation of lepton flavour universality, namely 

R{K) = ± 0.036, (1) 

in the range 1 GeV^ < < 6GeV^. This measurement 

is in tension with the theoretically clean SM prediction 
^SM(Ar) = 1.0003 ± 0.0001 [HI by 2.6cr. Combining 
these anomalies with all other observables for b —>■ sp'^p~ 
transitions, it is found that a scenario with NP in Cg^ 
(corresponding to the operator Ppbp'yvP) but not in 
Cg® is preferred compared to the SM by 4.3cr [TB] . 

Hints for lepton flavour universality violating NP also 
comes from the BaBar collaboration that performed an 
analysis of the semileptonic B decays B —>■ m- 

Recently, these decays have also been reanalyzed by 


BELLE |18| and LHCb measured B D*tv [19] . In 
summary, these experiments have found for the ratios 
B(bW) = Br(B ^ B)WT^)/Br(B ^ D<^*'>eiz): 


B(B)BaBar = 0.440 ± 0.058 ± 0.042 , (2) 

i?(^)BELLE = 0.375t°:“4 ± 0.026, (3) 

B(B*)BaBar = 0.332 ± 0.024 ± 0.018, (4) 

R{D*)belle = 0.293t°;°39± 0.015, (5) 

i?(A>*)LHCb = 0.336 ±0.027 ±0.030. (6) 


Here the first (second) errors are statistical (systematic). 
Combining these measurements one finds [20] 


R{D)exp = 0.388 ±0.047, 

R{D*)exp = 0.321 ±0.021. (7) 

Comparing these measurements to the SM predic¬ 
tions pT] 

Rsm{D) = 0.297 ±0.017, 

RsuiD*) = 0.252 ±0.003, (8) 


we see that there is a discrepancy of 1.8 cr for R{D) and 

3.3 cr for R(D*) and the combination corresponds approx¬ 
imately to a 3.8 cr deviation from the SM (compared to 

3.4 CT taking into account the BaBar results only [HD- 

Numerous models have been proposed in order to ex¬ 
plain the anomalies in & —>■ sp^p~ transitions (see for 
example Refs. [22U32] for Z' models and Refs. [33| [34] 
for models with leptoquarks) and the deviations from the 
SM predicitons in tauonic B decays [35H45] . 

Alternatively, a model independent approach using 
higher dimensional operators has been employed, as in 
the model independent hts [ilHlISI. In this context, it 
has been argued that as R{K) violates lepton flavour uni¬ 
versality (LEU) also lepton flavour could be violated in 
B decays [45] which might be linked to neutrino oscilla¬ 
tions [3B| . [55] While [IS] considered the effect of operators 
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at the B meson scale which are invariant under electro¬ 
magnetic gauge interactions only, also operators invariant 
under the full SM gauge group gTlIlH] have been consid¬ 
ered in Ref. . [53] Here it has been claimed than an 

simultaneous explanation of R{K), R{D) and R{D*) us¬ 
ing gauge invariant operators with left-handed fermions 
is possible [531153]. For this purpose, it was assumed that 
in the interaction eigenbasis only couplings to the third 
generation exist gsiEn] ( or are enhanced by com¬ 

pared to the second one [52]), while all other couplings 
are generated by the misalignment between the mass and 
the interaction basis (or are suppressed by small lepton 
mass ratios [53]). 

In this article we reconsider the possibility of explain¬ 
ing B —>• TV and the b —> fj, data with higher di¬ 

mensional gauge invariant operators, taking into account 
the constraints from B —>■ K^*^vv and using the results 
of the global fit to 6 —>■ s/r/r transitions. We extend the 
analysis of Ref. [53] and consider the possibility of lep¬ 
ton flavour violation (LFV) and compared to Ref. [35] we 
include the correlations due to SU (2) l gauge invariance 
and give quantitative predictions for B —>■ and 

Bs r/r. 

The outline is as follows: In the next section we col¬ 
lect the necessary formulae for the flavour observables. 
Sec. [Illj discusses the gauge invariant higher dimensional 
operators relevant for our analysis and Sec. jlVj presents 
our numerical results. Sec. fV] briefly reviews some possi¬ 
ble UV completions. Finally we conclude in Sec. jVIj 


II. FLAVOUR OBSERVABLES 
A. 6 —>■ sfj,^fj.~ transitions 


B. B^K^*'>vv 


Following Ref. m we write the relevant effective 
Hamiltonian as 

nT = (e’ot + c«o«) (12) 

= ^[sl^PL,Rb][v^l^, (1 - 7®) Vj ], (13) 

and « —1.47/s^. In the limit of vanishing right- 
handed sh current, the branching ratios normalized by 
the SM predictions read 


R 


KM 


1 ^ 

5 I: 


i,j = i 



(14) 


The current experimental limits are < 4.3 [53] and 
< 4.4 [55]. 


C. R D^*\v 

The effective Hamiltonian for semileptonic b M c tran¬ 
sitions is 

Heff = ^V,bC^,^[cj>^PLb][ia^PLiyj] , (15) 

with = Sij (for massless neutrinos) taking into 

account only left handed vector currents. In this case 
the ratios of branching ratios are 


b —>■ s£i£j transitions are defined via the effective 
Hamiltonian 

a=9,10 


o: 


9(10) 


a 

47r 




3\ ’ 


(9) 


where the primed operators are obtained by exchanging 
L Mr R. 

Concerning B M K* fj ,'^, B^ M and B M 

/B M Ae+e”, as already noted in Ref. [331153] . 
C'g'^ < 0 and = 0 is preferred by data. However, 
also the possibility Cg^ = —Ciq < 0 gives a good fit to 
data. Using the global fit of Ref. mm we see that at 
(1 ct) 2(7 level 


- 0.53(-0.81) > > (-1.32) - 1.54(10) 

-0.18(-0.35) > > (-0.71) - 0.91(.ll) 

Interestingly, the values of Cg^, C^g favoured by R{K) 
and B M K* lie approximately in the same 

range. [73] Furthermore, a good fit to the current data 
does not require C'g^^, hence in the following we neglect 
operators with right-handed quark currents for simplic¬ 
ity. 


R(Z)t*))EXP 

R{DM)sm 


Elc: 


i=i 


cb 

L 3j 


E 

i=i 


r^cb 

£ 


2 


(16) 


with i = e, II which has to be compared to Eq. ([^ and 

Eq. 0. 


D. Lepton-flavour violating B decays 

Here we give formulas for the branching ratios of LEV 
B decays following the analysis of Ref. [56]. We take 
into account only contributions from the operators 

and while neglecting contributions from operators 

with scalar currents not relevant for our analysis. For 
Be M (with £ ^ t') we use the results of Ref. [57] 
neglecting the mass of the lighter lepton. The branching 
ratios for B M K^^^ix^^B M K^'^ir^e^ are computed 
using form-factors obtained from lattice QCD in Ref. [55] 
(see also Refs. [HISl]). The final results read 
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Br [B, i+e'-] 

Bi[B -)■ 


tb^MbJI^ 2 


327r3 




= 10 


-9 


O-KMW 


cf + cf 


■ bK(-)U' 


Ml 

Olo -h Oio 


Cf - c. 


'Hi' 


2 2 
I r^ii' f-ilii' 

-h 0^0 ^10 


2 2 
~r CK*U' ^9 ~ ^9 -r(iK*i£' ^lo ~ ^10 


(17) 


with 


££' 

a-KU' 

bKW 

aK»u' 

bK*U' 

Ck*U' 

dK‘U' 

re 

lie 

9.6 ± TO 
15.4 ±3.1 

10.0 ± 1.3 
15.7T3.1 

3.0 ±0.8 
5.6 ± E9 

2.7±0.7 
5.6 ± E9 

16.4 ±2.1 
29.1 ±4.9 

15.4 ± 1.9 
29.1 ±4.9 


The formula for the branching ratio of Bg —>■ £~^f~ is 
symmetric with respect to the exchange of O 

, while in the case of il —>■ this sym¬ 

metry is broken by lepton-mass effects. There is a 
small difference between the theoretical prediction for 
the charged mode B^ —?> and the neutral 

one B^ —)■ due to the different B-meson life¬ 

time tb which we neglected fixing the numerical value of 
tb to the one of the neutral meson. Note that the results 
above are given for final states and not for the sum 

= £“£'+ + £^i'~ to which the experimental con¬ 
straints apply |60| . The only channel with r/i final states 
for which an experimental upper limit exists is 

Br [B+^ < 4.8 X 10"^ (18) 


neutrinos) we find for the terms relevant for the processes 
discussed in the last section (before EW symmetry break¬ 
ing) 


C {£il^PL£jdkl^tPLdi + n'y^^PLVjdk-f^PLdi) + 

^(3) 

—^ {2£i'y^PLVjUk^i,PLdi - i^a^PLVjdkj^iPLdi 

+ hrPLijdkl^PLdi) , (21) 


where are the dimensionless coefficients of the op¬ 

erators of Eq. (201. After EW symmetry breaking the 
following redefinitions of the fields are performed in or¬ 
der to render the mass matrices diagonal 


III. GAUGE INVARIANT OPERATORS 


di —>■ D^dL, ul —t U^ul, £l —>■ L^£l, v —>■ . (22) 


As we have previously seen, a scenario with left-handed 
currents only gives a good fit to data, cf. Eq. (11). In 
such a scenario SU{2)l relations are necessarily present. 
These relations are automatically taken into account once 
gauge invariant operators are considered. Therefore, let 
us focus on 4-fermion operators with left-handed quarks 
and leptons. There are two such 4-fermion operators in 
the effective Lagrangian 


TldimS = ^ OxCx , (19) 

where A is the scale of NP, which can contribute to 6 —>■ 
s££ transitions at tree-level imiiH]: 

= (L^L) (Qy^Q) , Qg) = (Ly'^uL) (Qy^r^Q) , 

( 20 ) 


We define for future convenience 


y(1.3)^(1.3)^1.3) ^ , (23) 


where are overall constants. Using constraints from 
the measured CKM matrix, i.e. V = WD, we finally 
obtain 


*-"9 ~ '-^10 


V2A^GFaVaV*g 

(jij __ ^ _ 

^ V2A^GFaVtbVA 


df + 


(AV.Y<;>fV> - 


A(3) a, 


(3) 

(^cb _ _^ ' 

^ V 2 A‘^Gf Vcb 


E 


(3) 

fc3 


^23 

^ij ^23 ) 1 


(24) 


where L is the lepton doublet and Q the quark doublet 
and the flavour indices are not explicitly shown here. 
Writing these operators in terms of their SU{2)b compo¬ 
nents (i.e. up-quarks, down-quarks, charged leptons and 


for the Wilson coefficients relevant for b sjjP ^ , 

B kP'i vv and B —>■ D^*\v respectively. Note that in 
the limit = C*^3) the contribution to i? —)■ 
vanishes. 
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FIG. 1: Allowed regions in the plane from b —>■ 

Sfj,^H~ data (blue) and B —> Ki/v (yellow) for ctsi, = 
ArcSin[I4i,] and A = 1 TeV. Note that here changing only 
has the effect of an overall scaling of A*-^A The contour lines 
denote Br[B —> K*Tfi] in units of 10“®. 


data simultaneously. From Eq. (27) we see that 


this ansatz is only possible for = C^^^but not if 
or (G(i) are separately different from zero. 

Therefore, we will focus in the following on scenarios 
with third generation couplings in the EW basis only, 
which correspond to a general rank 1 matrix in the mass 
eigenbasis, as suggested in Ref. In other words 

we have 


X = L^XL,Y = D^YD , X =Y = 


(30) 



Taking into account only rotations among the second and 
third generation one finds 


X = 

Y = 


0 

0 

0 

0 

0 

0 


0 

sin^ {a^r) 

-sin(a^.r)cos (a^r) 

0 

sin^ {usb) 

-sin(ast,)cos(Q;s&) 


0 

-sin(a^T-)cos (a^r) 


cos' 


(cT/tt) 


0 

-sin(ash)cos (a^fe) 


cos' 


(c^sb) 


Note that a rotation sin(as;,) ^ Vcb would require fine- 
tuning with the up sector in order to obtain the correct 
CKM matrix. 


IV. NUMERICAL ANALYSIS 

Since we have Cg^ = —C^q we find for the LEV B 
decays 

Br [B /Br [B « 1, (25) 

Br [B^ /Br [B K*t^h^] « 0.5 . (26) 

Therefore in the following, we will just present 
the numerical evaluation of Br [B —>■ while 

Br [Bg — 7 ^ and Br [B can be obtained 

by the appropriate rescaling. 

We also note that B ^ KvP imposes an upper limit on 
the absolute value of = —CIq and = —CIq valid 
for and separately. Neglecting the small NP 
contribution to and assuming no NP in the electron 
channel we find: 

V4.3x3/2«2.5, (27) 

< V3 X 4.3 X 3/2 - 2 + 1 « 5.2 . (28) 

This leads to the following upper limits valid in any 
model generating only or C'(i): 

Br[R ^ KT^J\ < 8.3 x 10"® . (29) 

However, this limit can be evaded for 
In Ref. 02] it was proposed that the MFV-like rela¬ 
tion Y 22 /A 33 = rri^/vn?^ could explain R{D^*'>) and b 


1. operator 

In this case we have neutral currents only. As a con¬ 
sequence, there is obviously no effect in R{D^*^), but 
b —>■ spX pb~ is directly correlated to H —>■ K^*^vv de¬ 
pending on the angle a^r- Note that a change in ash 
can be compensated by a change in and therefore 
does not affect the correlations among B K^*'>vv and 
b —> transitions. In Fig. the regions favoured 

by 6 —>■ spL^pL~ (blue) and allowed hj B —b KvD (yellow) 
are shown together with contour lines for B —> K*TpL in 
units of 10"®. Note that B Kvv rules out branching 
ratios for B —>■ K*Tpi above approximately 1 x 10"® and 
that the constraint from B —>■ Kvv^ being inclusive in 
the neutrino flavours, is independent of 

2 . Q\^ operator 

Here we have also charged currents that are related 
to the neutral current processes via CKM rotations. In 
Fig. i the regions allowed by H —)■ Kuiy (yellow) and 
giving a good fit to data for b —>• (blue) and (at 

the 2(7 level) for B —)■ D*tv (red) are shown for dif¬ 
ferent values of A*-®A Note that b —>■ spL^pT data can 
be explained simultaneously with R{D^*^) for negative 
0 (1) values of A^®^ without violating the bounds from 
B —>■ Kvv. Again, in the regions compatible with all ex¬ 
perimental constraints, the branching rations of LFV B 
decays to rpL final states can only be up to « 10"®. 
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FIG. 2: Allowed regions in the anT-otsb plane from B —>■ Ki^v (yellow), R{D*) (red) and b —^ (blue) for A = 1 TeV 

and = -0.5 (left plot), A^^' = -1 (middle) and A^®^ = -2 (right). Note that Usb = 7r/64 roughly corresponds to the angle 
needed to generate Vcb and that if A^^^ is positive, R{D*) and b —>■ cannot be explained simultaneously. 


3. and with A*-^^ = A^®^ 

In this case the phenomenology is then rather similar 
to the case of only. The major differences are that, 
as already mentioned before, the bounds from B ^ Kvv 
are evaded and the relative contribution to 6 —>■ s/i/i com¬ 
pared to R{D^*^) is a factor of 2 larger. In Fig. we 
show the analogous plot to the central panel of Fig. 
(A^^^ = = —1) for this scenario. Note that again 

R{D^*'>) rules out very large branching ratios for lepton 
flavour violating B decays in the regions compatible with 
b —data. We also consider the MFV-like ansatz 
[5^ with additional flavour rotations (light blue) which 
however differs only slightly for the ansatz with third 
generation couplings. 

V. UV COMPLETIONS 

Let us briefly discuss UV completions which can give 
the desired coupling structure. As discussed previously, 
the 4-Fermi operator is relevant both for R{K) and 
R{D^*'>). If is mediated by a single field, then there 
are only four possibilities: (i) Vector boson (VB) with 
the SM charges {SU{?>)c, SU{2)l,U(1)y) = (1,3,0), (ii) 
Scalar leptoquark (SLQ) with (3,3,—1/3), (iii) Vector 
leptoquark (VLQ) with (3,1,2/3), and (iv) Vector lepto¬ 
quark with (3,3,2/3). The vector boson (1,3,0) induces 
(S') 

only . On the other hand, the leptoquark fields result 

in particular combinations of and j52]. With 
the assumption of the third generation coupling, the rel¬ 
ative size of the effective couplings and the signs 

are determined as 

VB(1,3,0) : both positive and negative, (31) 


SLQ(3,3,-l/3) : A^^) = SA^^), A^^^ > 0, (32) 

VLQ(3,1,2/3) : A^^) = A^^^, A^^) < 0, (33) 

VLQ(3,3,2/3) : A^^^ = -3A(3\ A^^) > 0. (34) 

The coefficient Cq is proportional to A*-^^ -|- A*-^^ and a 

negative value is favoured by R{K). Therefore, the scalar 
leptoquark is rejected as a candidate. To explain R{D^*'>) 
simultaneously, A*-^^ itself must also be negative. This 
condition excludes the triplet vector leptoquark. If the 
experimental results are explained by the operator 
under the assumption of third generation coupling only, 
the possible mediators are the triplet vector boson or the 
singlet vector leptoquark. According to the analysis of 
the previous section, a good fit to flavour data requires 
a mediator mass of 0(1) TeV. This opens interesting 
prospects for the LHC, especially in the case of lepto- 
quarks that can be produced in proton-proton collisions 
via colour interactions and would decay to one lepton (r 
or more interestingly /r) and one jet (possibly a 6-jet). 


VI. CONCLUSIONS 

In this article we considered the effect of gauge in¬ 
variant dim-6 operators with left-handed fermions on 
6 —> B —> B —> B —>■ 

and Bs —> t^. For operators with left-handed quarks 
and leptons we find the correlations Br \B —)■ Kt^ « 

Br [B —>■ « 2Br [Bs —>■ We showed that 

the anomalies in 6 —>■ s/r/r data can be explained simul¬ 
taneously with R{D*). For this we considered scenarios 
in which third generation couplings in the EW basis are 
present only: A^^) ^ 0, A^^) ^ 0 and A^^) = A^^^ ^ 0. 
Taking into account A^^^ ^ 0 only, 6 —>■ data 

can be explained without violating bounds from B —>■ 
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FIG. 3: Allowed regions in the a^^T-asb plane from R{D*) 
(red) and b —>• sfi'^fj.~ (dark blue) for A = TeV and = 
= —1. The light blue region corresponds to the MFV-like 
ansatz for the lepton masses. Note that aab = 7r/64 roughly 
corresponds to the angle needed to generate Vcb and that the 
MFV-like Ansatz only differs marginally from the one with 
third generation couplings only in the region compatible with 
R{D). The contour lines denote Br[B —>■ K*Tfi] in units of 
10 “®. 


k(*) vv. However, in the allowed regions of parameter 
space, Br[i3 —)• can only be up to 1 x 10“®. In the 

case of A^^) 0, 6 —data can be explained simul¬ 

taneously with R{D*). In these regions Br[i? —>■ 
can again be only up to 10“®. Finally we considered 
A(3) = A(i) ^ 0. Such a scenario can be realized with a 
leptoquark in the singlet representation of SU{2)i^ (mak¬ 
ing an MFV-like ansatz for the lepton couplings possible) 
and constraints from B —>■ are avoided. Again, 

LFV B decays turn out to be of the same order as in the 
other scenarios. 


Note added — During the completion of this work, an 
article presenting a dynamical model with additional vec¬ 
tor bosons and third generation couplings appeared in 
which is generated m- 
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We denote the SM scalar particle predicted by Brout, 
Englert and Higgs as the ” Higgs particle”. 

Very recently, this discrepancy increased to 3.5 ct |62 |. 
Lepton flavour violating B decays in leptoquark models 
have been studied in and in Z' models in | 56 |. 

For an analogous analysis in the lepton sector see |631465| . 
Note that the fit to Eq. (Ill includes muon data only. 
However, as B > Ke'^e~ agrees rather well with the SM 
prediction, the effect on the global fit is expected to be 
small. Also the latest LHCb result for Bs —>■ 4>n'^ | 62| . 
which would slight increase the tension with the SM, is 
not included in the fit. 



